Quantum walks as a probe of structural anomalies in graphs 
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We study how quantum walks can be used to find structural anomalies in graphs via several 
examples. Two of our examples are based on star graphs, graphs with a single central vertex to 
which the other vertices, which we call external vertices, are connected by edges. In the basic star 
graph, these are the only edges. If we now connect a subset of the external vertices to form a 
complete subgraph, a quantum walk can be used to find these vertices with a quantum speedup. 
Thus, under some circumstances, a quantum walk can be used to locate where the connectivity of a 
network changes. We also look at the case of two stars connected at one of their external vertices. 
A quantum walk can find the vertex shared by both graphs, again with a quantum speedup. This 
provides an example of using a quantum walk in order to find where two networks are connected. 
Finally, we use a quantum walk on a complete bipartite graph to find an extra edge that destroys 
the bipartite nature of the graph. 

PACS numbers: 03.67.-a 



I. INTRODUCTION 

One of the most versatile quantum algorithms is the 
quantum search algorithm due to Lov Grover In its 
original form, it identified which Boolean function from 
a particular set was realized by a particular quantum or- 
acle. A Boolean function, f(x) maps n bit binary num- 
bers to either or 1, and the particular class of Boolean 
functions considered by the simplest form of the Grover 
algorithm are for all strings except one. We are given 
an oracle that realizes one of these functions; if we input 
x its output is f{x). Our task is to find which function it 
realizes, or, equivalently, for which input f(x) = 1, with 
as few calls to the oracle as possible. Classically one 
needs of order 2 ra calls whereas on a quantum computer, 
using the Grover algorithm, one needs only of order 2™/ 2 
calls. 

A variant of the Grover algorithm was defined for 
searches on graphs. First, one defines a quantum walk 
on a graph, which is a quantum version of a random 
walk 2] -[6]. Then the behavior of one of the vertices is 
changed so that it acts differently from all of the others. 
The object is then to find the distinguished vertex. This 
has been done for a number of highly symmetric graphs, 
such as the hypercube @, grids in different dimen- 
sions [H, [l(| , and the complete graph EH ■ The initial 
state of the walk cannot incorporate any knowledge of 
the distinguished vertex, and it is usually an equal su- 
perposition of all vertices, in the case of a coined walk, 
or an equal superposition of all edges, in the case of a 
scattering walk. The number of steps the walk must take 



in order to find the distinguished vertex is of the order 
of the square root of the number of vertices in the graph. 
Some of the latest studies of searches on graphs have fo- 
cussed on how the search is affected by the connectivity 
of the graph or by disorder in the graph [l2| , or searching 
in a graph in which there are several kinds of non-special 
vertices [l3j]. It should be noted that by constructing a 
quantum circuit that implements a quantum walk, these 
graph search problems can be rephrased as searches in- 
volving calls to an oracle. For an explicit example of this 
see [llj]. 

More recently, it has been found that quantum walks 
can find things besides distinguished vertices in a graph 
(l4j . In that study, walks on star graphs were examined. 
A star graph has a central vertex and N edges emanating 
from it, each of which is connected to its own vertex, so 
that the graph has a total of JV+1 vertices. We shall call 
the vertices besides the central vertex external vertices. 
If one adds an extra edge connecting two of the external 
vertices, it is possible to find the extra edge in approxi- 
mately y/N steps using a quantum walk. If one adds a 
loop to an external vertex, the result is the same, but if 
one adds a new vertex and an edge between that vertex 
and one of the external vertices, the quantum walk does 
not find the extra edge. So, it is unclear what kinds of 
structural anomalies can be found and what kinds can- 
not. 

Here we would like to continue our exploration of this 
subject. We begin by reviewing some of the results of 
14] and presenting more details. We then move on to 
several more examples. First we consider a star graph 
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with extra edges added connecting external vertices so 
that these external vertices form a complete graph. A 
complete subgraph of a graph is known as a clique. The 
idea is to use a quantum walk to find the vertices that 
comprise the clique. Note that what this does is allow 
us to find a part of the graph in which the structure 
of the network changes. In the star graph, the external 
vertices are only connected to each other through the 
central vertex. If we now form a region in which the 
density of connections increases, in particular in which 
the external vertices are directly connected to each other, 
we can use a quantum walk to find this region. Next, we 
consider two star graphs joined at one of their external 
vertices. In this case the walk starts on both of the stars, 
and we want to find the vertex where the two star graphs 
are connected. This shows that we can use a quantum 
walk to find where two networks are connected to each 
other. Finally, we look at the case of a complete bipartite 
graph. In this graph, the vertices are divided into two 
sets, and each vertex in one set is connected to all of 
the vertices in the other set by an edge, but no vertices 
within the sets are connected. Suppose we now add one 
edge connecting two vertices in one of the sets. We can 
use a quantum walk to find this edge faster than we could 
classically. 



II. STAR GRAPHS 



Throughout this paper we will be using the scattering 
quantum walk in which the particle "scatters" off the 
vertices of the graph [ll[ . There is another version of 
the discrete-time quantum walk, known as a coined walk 
0, 0]. In this type of walk, there is an extra system, 
the coin, that makes the step transformation unitary and 
controls the dynamics of the walk. The coined walk has 
been shown to be equivalent to the scattering walk, so 
that which one is used is a matter of preference [15| . 
We find the scattering walk more physically motivated, 
and it is the one we shall use. In this walk the particle 
making the walk sits on the edges of the graph instead 
of the vertices. Each edge has two orthogonal states. 
If the edge connects vertices j and fc, then one state is 
\j, k) corresponding to the particle going from j to fc, 
and the other is \k, j) corresponding to the particle going 
from k to j. The collection of all of these states, two 
for each edge, forms an orthonormal basis for the Hilbert 
space in which the walk takes place. In addition to the 
Hilbert space we need a unitary operator that advances 
the walk one step. In the scattering walk each vertex acts 
as a scattering center and is described by a local unitary 
operator that maps states entering the vertex to states 
leaving the vertex. The unitary operator that advances 
the walk one step, U , is simply made up of the action of 
all of the local unitary operators at the vertices. For a 
vertex, j, with n edges connected to it, we will generally 



use the operator 

n 

U\k,j) = -r\j,k) + t £ \j,l), (1) 

where r — (n — 2)/n and t = 2/n. This type of vertex 
behaves in the same way no matter from which edge it is 
entered. 

As was stated in the Introduction, a star graph has a 
central vertex, which we shall denote by 0, and N ex- 
ternal vertices, which we shall denote by 1 through N. 
Each of the external vertices is connected to the central 
vertex by a single edge. The dimension of the Hilbert 
space in which a walk on this graph takes place is 27V. 
The central vertex behaves as described in the preceding 
paragraph, and the behavior of the external vertices de- 
pends on the application. Since the walk on a star graph 
with an extra edge was discussed thoroughly in [TJ] , here 
we will describe what happens when we add loops to the 
external vertices. 

Let us first consider the case in which we add a loop 
to a single external vertex, which we shall take to be 
vertex 1. We shall denote the single state of the loop 
by |Zi). The operator U now acts on the states entering 
the external vertices as U\0, 1) = \l{), U\h) = 1 1 , 0) , and 
U\0,j) = \j, 0) for j > 2. Let us now define the states 
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and note that 



-r|0,l) 

|V>2) 



U\1,0) 
Ufa) 

U\ih) = t^W^ T|0,l)+r|^i). 



(3) 



Now, if S is the subspace spanned by the vectors 
{|0, 1), 1 1 , 0) , l^i), 1^2)}, we note that S is invariant 
under the action of U. This implies that if the initial 
state of the walk is in S, the entire walk will take place 
in S, which reduces the dimension of the space we have 
to consider from 2N to 5. This will be a feature of all 
of the problems we consider here, a drastic reduction in 
the size of the space due to the high symmetry of the 
graph. This type of dimensional reduction was first used 
by Krovi and Brun in studies of coined quantum walks 
[la ]. In this case, it means that U restricted to S, Us, is 
given by the 5 x 5matrix 
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where the basis is ordered as above in the definition of 
S. The state of the walk after n steps will be Ug \ipi n it), 
where |f/Wt) is the initial state of the walk and we shall 
assume that \ipi n it) <= S. In order to evaluate this, we 
want to find the eigenvalues and eigenstates of Us- 
The characteristic polynomial of Us is 

A 5_ rA 3 +rA 2_ 1 = ( A -l)(A 4 + A 3 +tA 2 + A + l) = 0. (5) 



we see that that the initial state is approximately equal 
to a superposition of two eigenvectors 

\TP m u) = ^(\v + )-\v-))+0(N-^). (10) 



If we now express —1 ± i\/t/3 = — exp(^fi6), where 
JtjZ we find that 



We see immediately that A = 1 is a root, but in order 
to find the others, we shall resort to perturbation the- 
ory. If N 3> 1, then t <C 1, and to find our zeroth 
order solution we set t equal to 0. The equation for the 
remaining roots then becomes (A + 1)(A 3 + 1) = 0, so 
that, to zeroth order, the remaining roots are —1 (twice) 
and exp(±i7r/3). Now we need to find the lowest order 
corrections to these eigenvalues. It turns out that the 
only interesting eigenvalue is — 1. This is because, as we 
shall see, the corrections to —1 are 0(N~ X / 2 ) while the 
corrections to the other eigenvalues are 0{N~ 1 ). In or- 
der to obtain a quantum speedup, we need the state to 
change substantially in 0(N 1 ' 2 ) steps. This will happen 
for a superposition of eigenstates whose eigenvalues are 
of the form Ao + O^ -1 / 2 ), where A is the zeroth or- 
der eigenvalue, but not for superpositions of states whose 
eigenvalues are of the form Ao + 0(N~ 1 ). In order to find 
the lowest order corrections to —1 we set A = — 1 + SX 
and substitute it back into the fourth order equation for 
A keeping only lowest order term. We find 



£HVw> 



SX = ±i 



3' 



(6) 



which is 0(N 1 / 2 ). The corresponding eigenvectors are 



with \v+) corresponding to —1 +i^/t/3 and \v-) corre- 
sponding to —1 — iy/t/'S. 

Now, for the initial state of our walk, let us choose the 
state 
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in — i 

-\/^ r (IV'i>-|V'2», 



which we can see is in S. Noting that 

IVi> - life) = i(K> - h>-». 



(8) 
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(11) 



Examining the form of U n \ipi n it), we see that when 
n6 = 7r/2 (this implies n is 0(N 1 / 2 )), the particle is 
either on the edge connected to the loop, with probability 
2/3, or on the loop itself, with probability 1/3. Now, in 
measuring where the particle is, we assume that we do 
not have access to the loop, otherwise we would know 
where it is, but we do have access to all of the edges. 
Therefore, if we measure the position of the particle after 
a number of steps satisfying n9 = 7r/2, we will with a 
probability of 2/3 find the particle on the edge connected 
to the loop. With a probability of 1/3, however, we will 
find no particle at all, and in that case we run the walk 
one more step, after which the particle will be on the edge 
connected to the loop. Therefore, by running the walk for 
0{N 1 / 2 ) steps, we have found which edge is connected 
to the loop with a probability close to 1. 

In comparing this procedure to a classical search for 
the loop, we shall assume that classically the graph is 
specified by an adjacency list, which is an efficient spec- 
ification for sparse graphs. For each vertex of the graph, 
one lists the vertices that are connected to it by an edge. 
This list can include the vertex itself, which means that 
there is a loop connected to that vertex. Searching this 
list classically would require O(N) steps to find the loop, 
while the quantum procedure will succeed in 0(y/N). 

The pattern of this calculation will be repeated for the 
other examples we discuss. First one finds an invari- 
ant subspace of small dimension in which the walk takes 
place. Next, one diagonalizes the unitary operator that 
advances the walk one step, U, in that subspace. This 
typically involves a perturbative approach to finding the 
eigenvalues and eigenstates. The zeroth order solutions 
are found by looking at the N — >• oo limit, and the small 
parameter in which one does the perturbation expansion 
is a power of 1/N. It is the eigenvalues that are degen- 
erate to zeroth order that lead to the interesting parts of 
the Hilbcrt space. One then identifies an appropriate ini- 
tial state, and calculates the action of U n on that state. 
Since this pattern holds for all of our calculations, we will 
present mainly the results in the body of the paper, and 
describe some of the details in the Appendix. 
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Before leaving the star graph proper, let us look at one 
more example, which was also discussed in [l4j|. Suppose 
that all of the external vertices except one, which we 
shall take to be vertex 1, have loops, and we want to find 
which vertex does not have a loop. Actually, we have 
to be a bit more careful in our description, because we 
are now going to assume we have access to the loops, 
so if there were one missing, we would know where it 
is. What we assume is that all of the external vertices 
are connected to loops, but the one connected to ver- 
tex 1 is a dummy loop. In particular, we assume that 
U\0,j) — \lj) and U\lj) = \j,0) for j > 2, and for vertex 
1, f7|0, 1) = e*>|l,0) and U\h) = \h). One only gets a 
quantum speedup for particular values of </>. 

This also reduces to a five-dimensional problem. The 
invariant subspace in this case is spanned by the vectors 
{|0, 1), |1,0), l^i}) IV^}}, where |^i> and \ip 2 ) are as 
before, and 



A' 



(12) 



We find that the characteristic equation of U restricted 
to the invariant subspace only has double roots in the 
]V4oo limit if <fi is tt, 7r/3 or — tt/3, and these are the 
values of <fi for which we obtain a quantum speedup. In 
the case that </> = it the appropriate initial state is given 
by 



N 



J2(\0,j) + \j,0) + \l j }), (13) 



and the particle becomes local ized on the edge with the 
dummy loop after n = {ir /2)y/z]i = O^N 1 / 2 ) steps. For 
4> = ±7r/3 different initial states are required, but the 
results are qualitatively the same. The details of the 
case 4> = t are given in the Appendix. 



III. STAR GRAPH WITH A CLIQUE 

Now suppose that we start with a star graph with N 
edges, and we add extra edges to it. The case of one 
extra edge was dealt with in [l4j], but now we wish to 
add enough edges so that a subset of the external ver- 
tices form a complete graph, or clique (see Fig. 1). In 
particular, we shall assume that vertices 1 through M 
form the clique, i.e. each of these vertices is connected to 
all of the other vertices in the set {1,2,... M} as well as 
to the central vertex. We will also assume that M -C N . 
This graph can be viewed as a network in which most of 
the participants are only connected through the central 
vertex, but there is a subset of participants who are di- 
rectly connected to each other. The object is to find the 
vertices in the clique. 




FIG. 1: A star graph with a clique, in this case a triangle. 
In this case, the operator U acts as follows 

N 

U\j,0) = -r\0,j)+t J2 l°< fc > 

k=l,k=£j 
M 

U\0,j) = -r\j,0)+i Y, \^ k ) forl<i<M 
U\0,j) = 1.7,0) forM+l<i<iV 

M 

U\j,k) = -r\k,j)+t\k,0)+t l fe >*> 

for 1 < j, k < M, (14) 

where r and t are as before, and r = (M — 2)/M and 
t = 2/M. We choose the initial state to be 



1 N 

|Vw) = ^=£(10,^-1^0)) 



(15) 



and after running the walk for 




2M- 1 
2M(M - 1) ' 



(16) 



steps, the particle is located on one of the edges connect- 
ing the clique and the central vertex with a probability 
of (2M - 2)/(2M - 1) (up to terms of order (M/N) 1 / 2 ) 
and it is located on one of the edges of the clique itself 
with a probability of 1/(2M— 1). We assume that we do 
not have access to the edges of the clique itself, so that 
when we measure the position of the particle we either 
find it on one of the edges emanating from the central 
vertex, or we don't find it at all, because it is on one of 
the edges of the clique. Note that the probability of the 
particle being on one of the edges of the clique decreases 
as the size of the clique increases. 

Classically one would have to search the adjacency list 
of the graph in order to find a vertex that is a member 
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FIG. 2: Two star graphs connected at one of their external 
vertices. 



of the clique, and one would have to check 0(M/N) ele- 
ments. This compares to the 0(y // M/N) steps the quan- 
tum walk must make in order to find one of the vertices 
in the clique. Once one finds one vertex in the clique, 
the rest are found by reading off the vertices adjacent to 
that vertex from the adjacency list in both the classical 
and quantum cases. 



IV. TWO STARS 

Now let us look at a different problem. We have two 
stars, each with N edges. They share one external vertex, 
so the stars are connected, but we do not know which one 
(see Fig. 2). The object is to find the shared vertex. 

Let us denote the central vertices of the two stars by 
A and B. In order to analyze a quantum walk on this 
graph, we shall assume that they share vertex 1 The ex- 
ternal vertices of the first star (with central vertex A) are 
{1,2, ...N} and those of the second star (with central 
vertex^) are {1, JV + 1, JV+2, .. .2N- 1}. The quantum 
walk in which we are interested takes place in an eight- 
dimensional invariant subspace spanned by the vectors 
|^) = \A,1), \i> 2 ) = \1,A), |i& 3 ) = \B,1), |V 4 ) = \1,B) 
and 



A' 



1^5) 
\<h) 



1 N 

3 — * 

1 2JV-1 

7W=i £ 

V j=N+l 
1 2JV-1 

= E ii.*>- 



y/N 



(17) 



These eight vectors form an orthonormal basis for the 
invariant subspace. All of the vertices behave as before 
except for vertex 1. We shall assume that the particle is 
transmitted at vertex 1, that is, there is no amplitude for 
it to be reflected there. That means that 



u\^x) = iv> 4 > um = 1^2). 



(18) 



The operator that advances the walk one step acts on the 
other basis vectors in the invariant subspace as 



Ufa) 

u\i> 6 ) 



-r|Vi) +WN- l|Vs> 
~r\^ 3 )+WN-l\^ 7 ) 



(19) 



where, as before, t = 2/N and r = (N — 2)/iV. 
We now start the particle in the state 



1 



2VN 



N 



^2(\A,j) + \j,A})-(\l,B} + \B,l)) 
3=1 



2N-1 

- (\B,j) + \j,B)) 

j=N+l 



(20) 



j=N+l 



that is, a superposition of all of the edge states in the first 
star minus a superposition of all of the edge states in the 
second. We find that after n = n^/N /2 steps, the particle 
is located with high probability (1 — 0(jV -1 / 2 )) on either 
the edge between vertices 1 and A or the edge between 
vertices 1 and B. Therefore, using the quantum walk 
we can find the external vertex the star graphs have in 
common with 0{sfN) steps, whereas classically we would 
have to search the adjacency lists of the external vertices 
of one of the stars, which means searching a combined list 
containing 7V+1 items (one item from each of the vertices 
connected only to the central vertex and two items from 
the vertex connected to both central vertices). Therefore, 
the quantum walk gives us a quadratic speedup. 



V. COMPLETE BIPARTITE GRAPH 

We will now consider a type of graph that is actually 
a generalization of a star graph. A bipartite graph is one 
in which the vertices are divided into two sets, and only 
vertices in different sets are connected by edges; there are 
no edges between vertices in the same set. A complete 
bipartite graph is one in which each element in one set 
is connected to all of the elements in the other set. In 
the case of a star graph, one set contains only the central 
vertex and the other contains the external vertices. We 
shall assume that there are Ni vertices in set 1 and N2 
vertices in set 2, so that there are N1N2 edges in total. 
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The vertices in set 1 will be labelled 1, 2, . . . TVi, and those 
in set 2 will be labelled iVi+1, + 2, ... N 1 +N 2 . Finally, 
we will add one more edge, between vertices 1 and 2. This 
edge destroys the bipartite character of the graph. What 
we want to determine is whether a quantum walk can 
help us find this edge. It could, in principle, be between 
any two vertices in set 1 or between any two vertices in 
set 2. However, we are going to analyze a situation in 
which the symmetry between the two sets is broken. In 
particular, we are going to assume Ni ^> N 2 , so that the 
extra edge is in the bigger set. So, in conducting a search 
what we are trying to do is to find an extra edge in set 
1. 

We need to define a quantum walk on this graph. 
There are now three sets of transmission and reflection 
coefficients. We have t x = 2/N 2 and n = (N 2 - 2)/N 2 
for the vertices {3, 4, . . . N{\, 

N x +N 2 

U\j,k) = -r 1 \k,j)+t 1 E l fc '*>, ( 21 ) 

i=JVi+l,i#j 

where Ni + 1 < j < Ni + N 2 and 3 < k < N u and we 
have t 2 = 2/Ni and r 2 = (N\ — 2)/N\ for the vertices in 
set 2 

iVi 

U\j,k) = -r 2 \k,j)+t 2 J2 IM>, (22) 
i=hift 

where 1 < j < JVi and Ni + 1 < k < N x +N 2 . Finally, we 
have the transmission and reflection coefficients for the 
vertices attached to the extra edge, t = 2/(N 2 + 1) and 
f = (N 2 - 1)/(N 2 + 1), 

JVi+JV 2 

U\j,l) = -f|l,j)+t|l,2)+t J2 l 1 ^)' 

;=jv 1 +i,; 7 y 
N1+N2 

U\j,2) = -f\2,j)+t\2,l)+t Yl l 2 '0, (23) 

J=JVi+l,J^j 

where N 1 + l<j<N 1 + N 2 and 

N!+N 2 

U\l,2) = -f|2,l)+t E l 2 ^) 

J=JVl + l 

N!+N 2 

U\2,l) = -f|l,2)+t E ILO- (24) 

J=JVl + l 

We find that in this case there is an invariant subspace 
of dimension five in which the walk takes place. 

Our next step is to find the characteristic equation 
for the resulting 5x5 matrix for U restricted to that 
subspace, and then, in order to find the zeroth order 
solution, take the limit as the number of vertices goes 
to infinity. Now, however, we have two parameters, Ni 
and N 2 , so there are different ways in which we could let 
the number of vertices go to infinity. As was mentioned 



earlier, we shall look at the case Ni — > 00 and N 2 fixed for 
our zeroth order solution and then calculate corrections 
to it. This result will correspond to the case Ni >• 1. 
We now start the walk in the state 



f 



V2^ViiV 2 



E E (\j,k)-\k,j)), (25) 

fe=l j=JVi + l 



and let it go for 



n=- y /N 1 (N 2 



2), 



(26) 



steps. We then find that the particle is on the extra edge, 
with a probability of N 2 /(N 2 + 2) and on one of the edges 
connected to the vertices linked by the extra edge with 
a probability of 2/(N 2 + 2). As usual, we assume that 
when we measure the position of the particle, we do not 
have access to the extra edge, so that if the particle is 
on the extra edge, we will simply not find it. Classically, 
the adjacency list for this graph contains 2N\N 2 + 2 el- 
ements, and since we know that the extra edge is in set 
1, we would only have to search half of them, i.e. the 
entries corresponding to the vertices in set 1. Quantum 
mechanically, after making approximately y/N\N 2 steps, 
our probability of ending up on an edge connected to one 
of the vertices connected to the extra edge is 2/(N 2 + 2), 
so in order to have a high probability of ending up on 
such an edge, we would have to repeat the walk approx- 
imately N 2 times, for a total number of steps of order 
N 2 y/NiN 2 . The ratio of the total number of steps in the 
quantum walks to the number of items in the adjacency 
list is approximately y/N 2 /Ni. Therefore, if Ni ^> N 2 
quantum walk gives us an advantage. 

Suppose, however, that we are faced with a different 
problem. We are given a complete bipartite graph that 
may, or may not, have an extra edge in set 1 that renders 
it no longer bipartite. What we would like to determine 
is whether there is an extra edge or not, and we don't 
care where it is. In that case, we only need to run the 
quantum walk a number of times of order one. If af- 
ter such a walk, we cannot find the particle, we know 
it is on the extra edge, so the graph does indeed pos- 
sess such an edge. This would require y/NiN 2 steps of a 
quantum walk. Classically we would still have to search 
the adjacency list, which has approximately NiN 2 items.. 
For this problem, the quantum quantum walk provides 
a greater speedup, the ratio of the number of steps of 
the quantum walk to the items in the adjacency list is 
(NiN^ 1 / 2 , but the classical procedure will give us more 
information. It will tell us where the edge is, while the 
quantum procedure will only tell us whether there is such 
an edge. 



VI. CONCLUSION 

We have studied a number of examples in which a 
quantum walk can find a structural anomaly in a graph 
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more efficiently than can a classical search. In most cases, 
the anomaly was an extra edge or a set of extra edges. 
The case of two stars was somewhat different in that 
there we found which two edges were linked, or, looking 
at it from the point of view of vertices, which vertex had 
two edges instead of one edge attached to it. 

These examples suggest that there is a class of graphs 
whose structure can be usefully and efficiently probed 
by quantum walks. What the general features of this 
class are, we do not know. It would be useful to identify 
features of a graph that would indicate that some of its 
properties can be ascertained by running a quantum walk 
on it. That remains a topic for future work. 

Recently there has been considerable experimental 
work on quantum walks on a number of different systems 
All but the last of these have been quantum 
walks on a line. The paper by Schreiber, et al. re- 
ported on an implementation of a two-dimensional quan- 
tum walk J22J . The rapid progress in this area leads us to 
hope that walks on more complicated geometries can be 
implemented, which would open the door to performing 
quantum walk searches. This could make some of the 
results presented in this paper accessible to experiment. 
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eigenstates are 



Appendix 



A. One missing loop 




for A = I + h/| 



for A = I — i\j |. 



We now choose 



and, setting = y/t/3, this yields 



U n \^ init ) = -= 



1 

3 



/ sin(n#) 
— sin(n#) 
v /273cos(n6») 
v/273cos(n60 
V y/2j3cos(n0) J 



0{N 



-l/2\ 



(29) 



1 N 

_ 53(10,^ + 1^0) + ^)) 

-^(K)-M) + o(tv- 1/2 ), (30) 



(31) 



From this equation we see that when n9 = ir/2, the par- 
ticle is, with probability close to one, located on the edge 
with the dummy loop. The cases cj> = ±7r/3 yield dif- 
ferent eigenvalues and eigenstates, and hence require dif- 
ferent initial states, but the results are qualitatively the 
same, the particle is with a probability close to one on the 
edge connected to the dummy loop after 0{N 1 / 2 ) steps. 



As was noted in the main body of the pa- 
per, this walk takes place in a five-dimensional in- 
variant subspace spanned by the orthonormal basis 
{|0, 1), j 1, 0) , \ipi), life)}- ln this basis, ordered as 
in the previous sentence, we have that 



B. Star graph with a clique 

This problem possesses a five-dimensional invariant 
subspace, S. Define 



U 8 = 



( -r ty/N- 1 \ 



WN - 1 r 

1 

\ 1 J 



The characteristic equation of Us is 



A s + re* A 3 + -rA 



i<t> 



(27) 



(28) 



which in the N — > oo limit becomes (A 3 — l)(A 2 +e 10 ) = 0. 
From this we see that 1 will be a double root of this 
equation if (f> = 7r, e 2 ^ 1 ^ will be a double root if <j> = 7r/3, 
and e~ 2m / 3 becomes a double root if <f> = — 7r/3. For the 
case (j) = 7r we find that the relevant eigenvalues and 



life) 

1^3) 

IV*) 
life) 



1 M 

M 



3=1 
M 



3=1 



M M 



y/M (M 



1 ^=T T \j>k) 



1 



y/N -M 

1 

y/N - M 



] = 1 k=l,k^j 
N 

j=M+l 
N 

j=M+l 



(32) 
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V 



WM - 1 

ty/M(N — M) 



These vectors are orthonormal and constitute a basis for 
S. With this ordering, the matrix for Us is 

tM - 1 ty / M(N — M) \ 

-f WM - 1 

r 

1 - tM 
1 / 

The characteristic polynomial for this matrix is 

A 5 + (t - 1)A 4 + [2(M - l)t + t- 2] A 3 
-[2(M- l)i + t- 2] A 2 - (f- 1)A - 1 =0, (34) 

which in the N — > oo limit becomes 

A 5 + (?-l)A 4 + (t-2)A 3 -(t-2)A 2 -(t-l)A-l = 0. (35) 

The N — > oo equation has a double root of —1. Setting 
A = — 1 + <5 A wc find that 



C. Two stars 

The dimension of this problem can be reduced still 
further, from 8 to 4. If wc define the vectors 



6X = ±J 2M S M ZlU±i6. 



(36) 



(2M - l)N 

The eigenvectors corresponding to these eigenvalues are 
/ 1 \ 



\v+) 



M 



2(2M - 1) 
for A = — 1 + i6, and 



-1/VM - 1 
- y(2M-l)/(2M-2 ) 
V iV(2M-l)/(2M-2) 7 



l«-> 



/ M-l 
2(2M - 1) 



-1/VM-l 
iV( 2M-l)/(2M-2) 
V -V(2M- l)/(2M-2) y 



for A 



= — 1 — iO. For the initial state we choose 
1 N 

IVw) = -=2(|o,j)-|j,o)) 



(38) 



= _(| u+ )-| t ;_))+0( v ^7iV). (39) 

We then find that up to terms of order (M/7V) 1 / 2 , 

/ 



U n \i> m u) = {-l) n 



M-l 



2M 



sin(n#) 
sin(n^) 
-( M-l)' 1 / 2 sm(nd) 

^)CQs(ng) 
^M^cosH) 



v -V 2 ^*^ 008 ^ y 

(40) 

Therefore, when n9 = it/2 we find that the particle is 
on one of the edges going from the central vertex to the 
clique with a probability of (2M - 2)/(2M - 1) and a 
probability of being on the clique itself of 1/(2M — 1). 



Wl) = 


T2 m 


- 1*» 


w 2 ) = 


>*) 


- 1«) 


w 3 ) = 


>*) 




w 4 ) = 


>*) 


- m) 



(41) 



then we find that the action of the unitary operator, U, 
is given by 



U\ Wl ) 
U\w 2 ) 
U\w 3 ) 
U\w 4 ) 



\w 3 ) 



-r\wi) + WN - l\w 2 ) 
r\w 2 )+WN- l\wi). 



(42) 



Therefore, if our initial state is in the subspace, S', which 
we define to be the linear span of the vectors = 
1, . . .4}, then the dynamics can be describe completely 
within this four-dimensional subspace. The initial state 



(37) IVw) = 



2y/N 



N 



J2((Aj) + \j,A))-(\l,B) + \B,l)) 

3=1 



2N-1 

- (\B,j) + \j,B)) 

j=N+l 



(43) 



can be expressed as 



I2N 



(K) + K)) 



N-l 
2N 



(\w 2 ) + \w 4 )), 



(44) 



so that it is, in fact, in S'. Therefore, we have reduced 
our problem to a four-dimensional one. 

We can go even further if we look at U 2 . Under 
the action of U 2 , the subspace S' splits into two two- 
dimensional subspaces, one spanned by \w 2 }} and 
the other spanned by {|ws), l^)}. In particular, we have 
that 



U 2 \w{) = r\wi) - WN - l\w 2 ) 
U 2 \w 2 ) = r\w 2 ) +WN -l\w!), 



(45) 



which means that in the {|u>i), \w 2 )} subspace U 2 can be 
described by the matrix 



U 2 



-ts/N 



_ ty/N- 

1 r 



1 



(46) 
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Similarly, we have that 



D. Complete bipartite graph 



Define the orthonormal set 



U 2 \w 3 ) = r\w 3 ) +WN - l\w 4 ) 
U 2 \w 4 ) = r\w 4 ) - WN - l\w 3 ), 



(47) 



which means that in the {IW3), |u>4)} subspace U 2 can be 
described by the matrix 



U 2 



WN - 1 r 



(48) 



So what we are left with are two two-dimensional prob- 
lems, both of which are, mathematically, equivalent to 
Grover searches. The eigenvalues of both matrices are 
A = r ± it\/N — 1, which we shall denote by e ±l6 , respec- 
tively. This implies that 9 ~ 2/^/N. It is straightforward 
to find the eigenvectors and to use them to raise the ma- 
trices to an arbitrary power. We find that 



U 2n \ Wl ) 
U 2n \w 2 ) 
U 2n \w 3 ) 
U 2n \w 4 ) 



cos(n9)\wi) — sm(n9)\w 2 ) 
sin(n9)\wi) + cos(n9)\w 2 ) 
cos(n9)\w 3 ) + sin(n6>)|w 4 ) 
— s\n(nff)\w 3 ) + cos(n9)\u>4}. 



(49) 



These equations imply that when n9 = n/2 our ini- 
tial state, \ipinit), will have been transformed into 
(l/\/2)(|wi) - | W3 )) up to terms of order A^ 1 / 2 . That 
means the probability that the particle is located on the 
edges where the stars are connected is almost one. The 
condition n9 = n/2 implies that n = n^/N/A, and the 
number of steps in the walk is twice that, or ny/N /2. 

Note that other initial states for this walk are possible. 
The state 



Winit) = ~r^\ Wl ) 



N-l 

N 



W2>, 



(50) 



which is an equal superposition of all of the outgoing 
states on the first star minus all the outgoing states on 
the second, will also lead to a successful search. After 
0(VN steps the particle will end up, to very good ap- 
proximation, in outgoing states on the connected edges. 



IV-i) 



^=(|1,2> + |2,1» 




Ni+iV 2 

E (b'^ + b'. 2 )) 

j=N!+l 

N1+N2 

E (iuhm) 

j=Ni+l 

JVi N x +N 2 

E E i;>*> 

3 j=N 1 + l 

E E im- 



V(iVi - 2)N 2 i 
1 

V(Ni-2)N 2 ^^ Ni:: 



(51) 



They span an invariant subspace, S, of U, the operator 
that advances the walk one step on a complete bipartite 
graph with an extra edge. The matrix of the restriction 
of U to S, Us is given by 



Us 



( -f ty/Ni \ 

-(r 2 -t 2 ) 2y/t^ 

t^JTh f 

2^F 2 r 2 -t 2 

V 1 / 



(52) 



The characteristic equation of this matrix is 

(A-l){A 4 +(f+l)A 3 +[f+l-t> 2 -t 2 )]A 2 + (f+l)A+l} = 0. 

(53) 

This problem has two parameters, N\ and N 2 , so there 
arc different ways to take the limit as the number of 
vertices goes to infinity. We shall consider the case N\ — > 
00 and N 2 fixed. This implies that to obtain our zeroth 
order solution we will let t 2 — > and r 2 — > 1. In this 
limit, the fourth order equation for A becomes 



A 4 + (f + 1)A 3 + 2f A 2 + (f + 1)A + 1 = 0. 



(54) 



We find that —1 is a double root of this equation, so 
we set A = — 1 + SX and substitute it into the actual 
characteristic equation keeping only the smallest terms. 
This gives us S A = ±i9 where 
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2t 2 



N 2 + 2 

The corresponding eigenvectors are 



(55) 



\v+) = 



xJ1{N 2 + 2) 



-2 y/(^ 2 + 2)/2 

V iy/(N 2 + 2)/2 J 



(56) 
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for A = — 1 + i9, and 



After n steps, the state of the system is 



\v-) 



^2{N 2 + 2) 



i^ {N 2 + 2)/2 
V - V(^2 + 2)/2 y 



(57) 



for A = —1 — Both of these expressions are valid up 

— 1/2 

to corrections of order N 1 

For our initial state we choose 



Ni N!+N 2 



|Vw) 



- 1 V V (\j,k)-\k,j)) 



V2 



K)-M) + o(7vr 1/2 ). 



(58) 



f/"IVw> = 



1^2 Sin ( n61 ' 

sin(n6>) 



■■ sin(n0) 



-^cos(n0) 



1 . cos(n9) 



(59) 



Note that when n6 = -k/2 the particle is on the extra 
edge, with a probability of iV 2 /(iV 2 + 2) and on one of 
the edges connected to the vertices linked by the extra 
edge with a probability of 2/(N 2 + 2). 
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